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The irradiation of a dilute cloud of cold atoms with a coherent light field produces a random 
intensity distribution known as laser speckle. Its statistical fluctuations contain information about 
the mesoscopic scattering processes at work inside the disordered medium. Following up on earlier 
work by Assaf and Akkermans [Phys. Rev. Lett. 98, 083601 (2007)], we analyze how static speckle 
intensity correlations are affected by an internal Zeeman degeneracy of the scattering atoms. It is 
proven on general grounds that the speckle correlations cannot exceed the standard Rayleigh law. 
On the contrary, because which-path information is stored in the internal atomic states, the intensity 
correlations suffer from strong decoherence and become exponentially small in the diffusive regime 
applicable to an optically thick cloud. 

PACS numbers: 05.60.Gg,42.25.Dd,42.50.Ct 


I. INTRODUCTION 

In typical light-scattering experiments, a monochro¬ 
matic light beam irradiates a medium with a certain 
wave vector and polarization (incident channel a), and 
the transmitted light intensity T ab is then measured for a 
certain final direction and polarization (scattered chan¬ 
nel b). When the sample consists of randomly located 
scatterers, one observes a complex intensity pattern with 
large fluctuations known as a speckle pattern [I] . It orig¬ 
inates from the interference between a large number of 
field amplitudes scattered by the sample. This pattern 
depends intricately on the positions of the scatterers, on 
the frequency and polarization of the incoming light, on 
the direction of observation, etc. 

Since light scattering experiments allow for a precise 
control of the incident beam and accurate analysis of 
the scattered beam, they have been successfully used for 
thorough investigations of theoretical concepts originat¬ 
ing in the field of mesoscopic electron transport, like weak 
and strong localization, and universal conductance fluc¬ 
tuations |2]. Already well documented for classical scat¬ 
terers, light-scattering experiments investigating coher¬ 
ent multiple scattering have entered a new realm when 
cold atoms were used as scatterers [3]. It was rapidly 
realized that one crucial difference between classical and 
atomic scatterers resides in the atomic quantum internal 
structure, namely, the Zeeman degeneracy of the atomic 
dipole-transition levels ® E]. Each atom behaves like 
a freely orientable magnetic impurity that can exchange 
angular momentum with the scattered photon and thus 
store which-path information [B]. Tracing out the atomic 
magnetic degrees of freedom then generally results in 
strong decoherence, as observed by a drastic reduction 
of the coherent backscattering signal |7j compared to the 


case of non-degenerate point scatterers [S] . 

To capture the generic features of a speckle pattern, 
one usually performs an average over spatial config¬ 
urations of the scatterers. Typical quantities of in¬ 
terest are then the average transmission T ab and the 
intensity-intensity correlation T ab T a i b i where the over¬ 
bar denotes the spatial configuration average. When¬ 
ever the transmission T ab originates from the coherent 
superposition of a large number of uncorrelated, com¬ 
plex random amplitudes, the central limit theorem as¬ 
serts that its probability distribution is the Rayleigh law, 
P(T ab ) = T~ b exp {-T ab /T ab } [2^ As as consequence, 
the statistical fluctuations obey = n\T ab . 

Does the Rayleigh law hold also for the intensity- 
intensity correlations of light scattered by cold atoms? 
This question was first addressed in Ref. [3], where it was 
predicted that the Zeeman degeneracy leads to a drastic 
increase of fluctuations in excess of the Rayleigh limit 
T 2 fe = 2 T ab . This claim raised a controversy [TUI TT1 . 
With this article, we present a detailed theoretical re¬ 
analysis of the protocol proposed in (9,. In short, we 
ascertain that the internal Zeeman degeneracy reduces 
short-range speckle intensity correlations instead of en¬ 
hancing them. 

The rest of the paper is structured as follows: In the 
next Sec. [TTJ we state the relevant photon-scattering am¬ 
plitude, determine the ensemble-averaged intensity, and 
derive the appropriate expression for the intensity cor¬ 
relations. In Sec. m we compute the diffusive intensity 
correlations analytically as a function of Zeeman degener¬ 
acy, finding that they decay exponentially with the sam¬ 
ple thickness for atoms with a degenerate internal ground 
state. Sec. m concludes. 
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II. SPECKLE CORRELATIONS AND THE 
RAYLEIGH CONSTRAINT 

This section introduces the transition amplitude for 
the scattering of a photon from immobile atoms with 
internal structure, and determines the intensity and its 
correlations upon ensemble averaging over internal and 
external configurations. It is shown that the interfer¬ 
ence of amplitudes visiting atoms on different scattering 
paths can only involve Rayleigh transitions that preserve 
the internal state. Conversely, elastic Raman scattering 
events act as ideal which-path markers and thus reduce 
the relative weight of the interference terms. Finally, we 
prove for the proposed setting that the relative fluctua¬ 
tions are rigorously bounded from above by the Rayleigh 
limit. 


A. Light-scattering amplitudes 

Consider a dilute cloud of N identical atoms, cooled to 
a temperature low enough such that they can be consid¬ 
ered immobile at fixed classical positions R = {r i\j = 
1 neglecting recoil effects and Doppler shifts. 

The density is low enough for quantum denegeracy effects 
to be negligible. Then, each atom is an ideal classical 
point scatterer for electromagnetic radiation. We look at 
light being scattered by a closed optical dipole transition 
between two internal electronic levels, with total angular 
momentum J g in the ground state and J e in the excited 
state. In the absence of external fields, these levels are 
degenerate, and the ground state is spanned by magnetic 
quantum numbers m in the range — J g < m < J g . The 
set of magnetic quantum numbers for all atoms is noted 
M = {mP\j = 1 We assume that the atomic 

sample is uncorrelated, and statistically homogeneous, 
in the following sense: 

(?) the spatial ensemble average (• ) = /ILdr7(R) 
over the classical positions is taken with the sepa¬ 
rable probability density 

/(R)=n/( ri ) 

i 

where generically /(r* ) = 1/V with V the sample 
volume accessible to each atom, and 

( ii ) the ensemble average (•) = tr{p-} over the magnetic 
quantum numbers uses the direct product state 

N 

P = (g)p\ (2) 

i=l 

where p l = X) TO » Pm i |m’)(m® | is the initial state of 
atom ?. 

In the present section, the initial populations 0 < p m i < 1 
are left arbitrary. Only for the analytical calculations in 


the following Sec. Ill all atoms are assumed to be in the 
completely isotropic state with p m i = 1/(2 J g + 1). 

Consider now the scattering of a monochromatic light 
beam by the atomic sample. An incident field mode with 
wave vector k a and transverse polarisation e a is scat¬ 
tered into a superposition of final modes kf,,ef,. The 
incident light intensity is small compared to the satu¬ 
ration intensity of the atomic transition, such that only 
the linear response, i.e., scattering of a single photon at 
a time, has to be considered. As a consequence, the light 
is scattered elastically, |k a | = |k{,| = k. The internal 
state of an individual atom, however, couples to the po¬ 
larization of the scattered photon and thus can either 
remain the same (Rayleigh transition) or change (degen¬ 
erate Raman transition), with the concomitant change in 
the photon polarization imposed by angular-momentum 
conservation. The scattering amplitude from the initial 
state |k a , e a , M a ) to a final state |k{,, eb, Mb) is the matrix 
element 


A ab = (M b \T ab \M 0 ) (3) 

of the transition operator (also called T-matrix [12]) 
T a b = (kh,Ch|T(R)|k a ,e a ) acting on the internal states 
of all atoms [13! . 

Throughout the paper we assume that the atomic 
cloud satisfies the diluteness conditions n\ 3 <C 1 and 
fcl > 1, where n is the number density of atoms in the 
cloud, A = 2-7T /k the light wavelength, and l the scatter¬ 
ing mean free path. In this case the light propagation 
can be described semiclassically in terms of scattering 
paths I = {i\,i 2 ,... }, defined as the ordered set of light- 
scattering atoms at positions {r*|? £ 1} =: R J . The 
probability amplitude © then is the coherent superpo¬ 
sition 

A a b = A^b- ( 4 ) 


The diluteness condition also implies that one may safely 
discard recurrent scattering events, in which scatterers 
are visited repeatedly. Therefore, each atom appears at 
most once in any given path I. 

Importantly, only an atom ? £ / in the actual scat¬ 
tering sequence can have its magnetic quantum number 
to* affected. All other atoms j £ J c = {j\j ^ 1} are 
bound to be mere spectators in this scattering sequence, 
and their magnetic quantum numbers remain unchanged. 
In mathematical terms, this is expressed by writing the 
transition operator for path / in the form T 1 (§>1/ , where 
T 1 = T(R J ) operates only on the scattering atoms, and 
1 J does nothing to the atoms in the spectator ensem¬ 
ble. By partitioning the magnetic quantum numbers 
M = M 1 U M 1 over scatterers and spectators, we can 
write the scattering amplitude for the path / as 


Kb = (M^tiblMi) (Ml 


M, 


i c 


(5) 


The second factor imposes the constraint to/ = to/ for all 
spectator atoms j £ J c . An atom whose quantum num¬ 
ber is not changed during the sequence I can either be a 
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member of 7 and make a Rayleigh transition ( rri * l b = m l a ), 
or it can be a member of the spectator set 7 C . But it is 
immediately clear that any atom that makes a Raman 
transition (m| 7 ^ m l a ) must belong to a scattering se¬ 
quence 7. In other words, Raman transition events are 
path markers. 


B. Intensity 

The transition probability for the scattering process 
a — > b is \A a b\ 2 = A a bA^ b . Each of the amplitudes A ab 
here is the superposition Q of path amplitudes, such 
that the total probability splits into a sum of diagonal 
(7 = J) and off-diagonal (7 7 ^ J) contributions: 

\A a b\ 2 = E AibAil + E AibA'ab- (6) 

1 ijtj 

Because of the interference terms, this probability will in 
general depend quite sensitively on the atomic positions, 
and on the light held vectors (k a , e a , k^, e&); the resulting 
random light intensity distribution is known as an optical 


speckle pattern. In the following, we first discuss the 
effect of the internal degrees of freedom, for a fixed spatial 
configuration, and then proceed to spatial averages. 

1 . Internal ensemble average 

The path amplitude © depends on the internal con¬ 
figurations, M a and Mb . We assume that these internal 
quantum numbers are not under microscopic control and 
trace them out by averaging over the initial M a with 
weights pm 0 = Yl i Pm i and summing over the final Mb, 

T ab = Y PMjA ab \ 2 = (A ab A* ab ). (7) 

M a ,M b 

Inserting the double path sum ([h]) one thus arrives at the 
transition probability 

T ab = E (AgbAgl) + E (AlbA^h) . (8) 

1 i^j 

Taking into account the product structure ([5]), the inter¬ 
nal averages are 


(-A-lbA, 


ab / 


Y ( M « \ J ab\ M b ) ( M b 11 it | Mi ) 

M a ,M b 



Mf 



(9) 


There is an important difference between the diago¬ 
nal, or incoherent contributions (J = J) and the coher¬ 
ent contributions (7 7 ^ J) to the total intensity ([ 8 ]). This 
difference is best appreciated in the single-scattering ap¬ 
proximation (valid for an optically thin medium, such as 
a collection of a few scatterers), where the scattering path 
1 1 —> i reduces to the visited atom M- In this situation, 
using © with I = J = i, one sees that the incoherent 
contribution of atom i, 

( At ab A ab)= Y P m i ( TO “ I ^06 I m b) ( m b \ tab | m a ) (10) 

m a’ m b 

involves all transitions between its internal states that 
are allowed by the dipole selection rules. Notably, these 
can include degenerate Raman transitions (m b 7 ^ m[ J. 
The internal quantum numbers of the N — 1 spectator 
atoms are irrelevant and disappear by the elementary 
normalization P m i l( m al TO b)| 2 = tr {fP} = 1. In 

contrast, in the coherent contribution from two different 
atoms, 

(• KbKb ) = E P ^K( TO al^al| TO o)( m a|Ta6|m l 0 ), (11) 
m z a ,m 3 a 

each of the two atoms i 7 ^ j belongs to the spectator set of 
the other (i £ J c and j £ 7 C ). Consequently, two overlap 
factors (m J b \rn J a ) enforce Rayleigh transitions for 


both atoms i and j in the step from Eq. ([ 9 ]) to Eq. (Ill, 
and the N — 2 remaining atoms in their joint spectator 
set (I U J) c are traced out as irrelevant. 

The preceding reasoning, together with (101 and © 
written again with i I and j 1—7 J, holds also for 
multiple scattering sequences, as long as the two paths 

I, J do not share a common scatterer, In J = 0, such 
that I C J c and J C 7 C . This is actually the dominant 
situation in a dilute cloud with N 3> 1 atoms. For future 
use, we note that the coherent internal ensemble average 
then decouples as 


(A 


1 aJ*\ — 

b A ab) ~ 


T / 

1 ab 


J j 


(7n J = 0), (12) 


where (Tab) = t r {pT^ b } involves only atoms of path 7. 

It is a crucial insight for this paper that the internal 
path markers impose the Rayleigh-scattering constraint 
for the interference terms, whereas the incoherent terms 
now include also the Raman scattering. This obviously 
complies with the venerable rule of scattering theory (in¬ 
herited from the superposition principle of quantum me¬ 
chanics) that amplitudes for indistinguishable processes 
should be added coherently, whereas the probabilities 
for distinguishable processes must be added incoherently. 
And indeed, as already indicated above, as soon as a sin¬ 
gle atom makes a Raman transition, this event acts as 
a path marker that renders the amplitude distinguish- 
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able from all the others where this atom was not visited. 
Thus, the incoherent processes, including the Raman 
scattering events, are much more numerous compared to 
the non-degenerate case. As a consequence, the Rayleigh- 
scattering constraint reduces the apparent speckle con¬ 
trast with respect to this enhanced background. 

Of course, the reduction of interference visibility in the 
presence of which-patlr information is not a new insight 
special to mesoscopic systems, but has on the contrary 
been very well studied in numerous contexts. Promi¬ 
nent examples include the quantum-optical realization of 
Young’s double-slit experiment using two trapped ions 
with internal structure Hi fi5] . Mach-Zehnder photon 
interferometry with path encoding by polarization HU, 
and matter-wave interferometry with path encoding in 
internal electronic states [18; 'T9l. 


2. Spatial ensemble average 


Consider again the path-resolved transition probabil¬ 
ity © for light scattering in transmission across a di¬ 
lute cloud of scatterers, as depicted in the upper half of 
Fig. 1 [a). In the interference terms, all path pairs I ^ J 
involve large, random phase differences. (In a trans¬ 
mission geometry, one can neglect the pairing of time- 
reversed paths that is responsible for coherent backscat- 
ting. As pointed out before, one can also discard subdom¬ 
inant contributions where / and J share common scat¬ 
terers.) Under the spatial ensemble average over random 
atomic positions, the interference terms I ^ J vanish, 
and thus the average intensity contains only the paired 
paths / = J of the incoherent contribution, 


Starting from the connected intensity correlator 


5T ab 5T a i b i = T ab T a f b f - T ab T a ' b >, 


(14) 


one defines the magnitude of fluctuations compared to 
the background, 


^ 5T ab 5T a , b , 

( -'ab,a'b' — — — 

J- ab-l a'b' 


(15) 


as 


^ the quantity of central interest. The product of inten¬ 
sities 0 before spatial averaging reads 

TabTa'b'= Y ( A lbA J J) (A« bl A%) . (16) 

IJKL 

Along with Ref. [5], we assume here that the atoms are 
in the same, completely uncorrelated state ([2]) at the ar¬ 
rival of each pulse, such that the averages over the inter¬ 
nal configurations decouple. Fig. [lja) represents these 
amplitudes schematically. Just as for the average inten¬ 
sity (13), the only terms surviving the spatial average 


are those where the phases accumulated between differ¬ 
ent atoms are prefectly compensated. This leaves two 
possible path pairings, namely (/ = J,K = L ) and 
(J = L,J = K), depicted in Fig. 0b) and (c), respec¬ 
tively. 

Following [SJ, we assume that the paths / and K do 
not share any scatterer (In AT = 0).[2O; The first pairing 
(/ = J,K = L) then yields the product of the average 
intensities, 


E (KbAab) E ( A a'b' A % ) = TabTa'b'- (17) 


K 


As a consequence, the speckle correlations (141 are en¬ 
tirely due to the second pairing (J = L, J = K)\ 


Tab = Y{ V l V ab)’ ( 13 ) 

I 

as depicted in the upper half of Fig. ljb). This average 
transmitted intensity in the far field is a rather struc¬ 
tureless, smooth function of wave vectors and polariza¬ 
tion vectors without mesoscopic signatures of coherence. 
More interesting information can be gained from inten¬ 
sity correlations, to be discussed next. 


C. Intensity correlations 


We now consider the protocol proposed in [9] and an¬ 
alyze the correlations between the two speckle patterns 
obtained in consecutive scattering sequences a —> b and 
a' —> b ', recorded one shortly after the other. The wait¬ 
ing time between pulses is supposed to be short enough 
for the atomic positions to remain unchanged. Then, 
after each double pulse, the atoms slowly move to new 
positions. This two-pulse experiment is performed many 
times, until the ensemble average (•) over random atomic 
positions is achieved. 


ST ab ST a>b> = Y ( Al ab A«*) (A^A%). (18) 


IK 


Since I fl K = 0, both external and internal averages 
decouple. In particular, Eq. (12) can be used, resulting 
in 


sT ab 6T a/b/ =Y( ty(+£,*) EE^'Kw ( 19 ) 

i 

Yi^iU'b) ■ ( 20 ) 


The derivation of this expression is the key achievement 
of our paper. In essence, we find that under the as¬ 
sumption of uncorrelated distributions, the mesoscopic 
speckle correlations are given by a path-averaged prop¬ 
agator involving products of independent interna l aver¬ 


ages (Tf lb )(T^, fc ,)*. As explained in Sec. IIB1 above, 


this is ultimately a consequence of the Rayleigh scat¬ 
tering constraint for interference contributions. Ref. 0 
predicted unphysically large speckle correlations essen¬ 
tially because the Rayleigh constraint was not taken into 
account; herein lies the difference between Eq. (5) in the 
Reply by Assaf and Akkermans HU and our result (|20[). 
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FIG. 1. (Color online) Amplitudes for light scattering from mode a = (k a , e a ) to b = (kb, e b ) in a first pulse and from a' to b' 
in a second pulse, across a dilute cloud of immobile point scatterers. (a) Possible paths I, J, K, L, drawn with 4 representative 
scatterers each, (b) Pairing I = J and J = K, yielding the product T ab T a /y of average intensities, (c) Pairing I = L and 
J = K, describing intensity correlations 5T ab 8T a / b i . 


D. The Rayleigh bound 

We proceed to prove on general grounds that under the 
stated hypotheses (disjoint paths I n K = 0 and uncor¬ 
related, but otherwise arbitrary distributions of internal 
and external degrees of freedom), the relative fluctuations 
(15) are bounded from above by the Rayleigh law, 


Cab,a'b' — 1 - 

Proof: In a first step, let 


E t 


T!.„ = 


-fi.t'), 


( 21 ) 


( 22 ) 


denote a scalar product with respect to the trace over 
the atomic external degrees of freedom. The speckle cor¬ 
relations (20) are then bounded by the Cauchy-Schwartz 
inequality 


5T ab 5T a , b , = 


(t,t')| 2 < (t,t) (t',T'). (23) 


Equivalently, the general correlator (151 cannot be larger 


than the geometric average of the diagonal correlators 

Cab,ab = \{T,t)/T ab \ 2 : 


C a b,a'b' — \JC abi a b (J a ' b i , a ' b '. 


(24) 


In a second step, we prove that each diagonal correlator 
obeys the Rayleigh bound, 


C a b ab — 


Ei {jib) 


T I 

ab 


EAUlZb) 


< l. 


(25) 


This bound is a corollary of the Cauchy-Schwarz inequal¬ 
ity 

(26) 


now for the scalar product {A^ B) = tr {pA^B} of the 
trace over the internal degrees of freedom. Setting A = 1 


and B = T, one has 


T 1 

1 ab 


< 


j/f-r/ 

ab ab 


(27) 


and by linearity of the path sum and ensemble average, 


E 


fi 

ab 


<E( f a l^b 


(28) 


whence (251. Together with (|24j) , this yields the Rayleigh 
bound ( |21[ ). 


Clearly, expression (251 shows that a necessary and suf¬ 


ficient condition for satisfying the Rayleigh law C a b,ab = 
1 is that the internal ensemble average factorizes, 
<t'lT' 6 > = |<t^)| 2 . And we know from Ref. [5]— at 
least for the isotropic internal state considered in [5] ■ 
that the average intensity is equal to the averaged am¬ 
plitudes squared if and only if the ground state is non¬ 
degenerate (J g = 0). Therefore, we can formally con¬ 
clude that the fluctuations reach their maximum value, 
C a b,ab = 1 , only for J g = 0. In all other cases J g > 0, the 
fluctuations are reduced, C a b,ab < 1 - 


III. DIFFUSIVE INTENSITY AND SPECKLE 
CORRELATIONS 

In the following, we re-calculate the speckle correla¬ 
tions along the lines of [5] , but with the ensemble average 
over internal quantum numbers as defined in Eq. (201. 


We first recall how to compute the average light inten¬ 
sity scattered by a single atom with internal degeneracy 
in the isotropic ground state 


1 


Po = 


2J„ 


1 ^ 


% 


\{A^B)\ 2 < (A^A){B^B), 


(29) 
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Then we determine the single-scattering vertex for the 
correlation propagator. In a second step, we compute 
the diffusive intensities and their correlations. Finally, 
we determine the relative fluctuations observed in trans¬ 
mission across an optically thick slab. 


A. Single-scattering vertices 

The dipole scattering of a photon from mode (k a , e a ) to 
mode (kb, e b ) by a single atom at position r* is described 
by the scattering operator 


the single-scattering intensity vertex (331 which couples 


components k = 0,1, 2—and results in 

= M eg (e a ■ e* b )(e b , • <,)■ (35) 

The scalar products between polarization vectors here 
mean that the atom radiates exactly the polarization it 
has absorbed, in both sequences a b and a' >->• b', as 
required by angular-momentum conservation under the 
Rayleigh-scattering constraint. In this sense, the speckle 
correlations described by Eq. (201 originate from a much 


tab = ( e 6 ' d)(d ■ e a )e i ( ka-k| ’- , ' r< . 


(30) 


more restricted set of events than the full intensity. We 
are thus entitled to anticipate that internal degeneracies 
will rather reduce the observable speckle correlations, in¬ 
stead of enhancing them. 


Proportionality factors that are irrelevant in the present 
context have been set to unity, and d is the reduced 
atomic dipole operator connecting the ground-state man¬ 
ifold with angular momentum J g to the excited-state 
manifold with angular momentum J e . For the isotropic 
internal state (29), the average single-scattering inten¬ 
sity tr{p(,t^t^} = (t^t J a6 ) can be evaluated in closed 
form using the techniques of irreducible tensor operators 
0 Eq. (39)]: 

= M eg [uq|e a • el\ 2 + W 2 \£ a • 6b|“ + w 3] ■ 


Here, 


M = 2J e + ! 
eg 3(2 Jg + 1) 


(31) 


(32) 


denotes a ratio of level degeneracies, and the weights Wj 
are rotational invariants that depend only on the total 
angular momenta J g , J e . It is advantageous to define the 
reduced vertex function I(e a ,e* b ,e bl e* a ) = 
given by 

H e a, e b> e b, c*) = Wi\e a ■ e * b \ 2 + w 2 \e a ■ e b \ 2 + w 3 . (33) 

This function is essentially the radiation pattern, propor¬ 
tional to the differential scattering cross section, in the 
given polarization channel e a —> e b . 

Let us now examine the single-scattering building 
block for the correlation propagator (201, namely 


^U= M eg 1 ( T l ab )(TU 


(34) 


In the bulk of the scattering medium, one has k a = k a ' 
and fo, = kb' [cf. Fig. [ljc); the dependence on the exter¬ 
nal momenta will be re-established below], and we can 
focus on the dependence on polarization vectors. Be¬ 
cause of the Rayleigh-scattering constraint, the correla¬ 
tion vertex (341 is proportional to the product of two 
independent traces over internal states. Each trace se¬ 
lects the scalar component k = 0 of the dyadic tensor 
operator dd —a considerable simplification compared to 


B. Multiple scattering: the diffuson 

In order to sum the multiple-scattering intensity se¬ 
ries, it has proven useful m to write the cartesian ele¬ 
ments Iu-jk = w 1 5 i j5 k i + W 2 Sik 8 ji + w 3 SuSjk of the single¬ 
scattering 4-point vertex (33) as a sum of its irreducible 


tensorial components in the incident and scattered po¬ 
larization channels, respectively: 


lil-jk — £ ^nTu-jk- 


(36) 


K— 0 


The projectors onto the three subspaces k = 0,1,2 of 
dimension 2k +1 describing population, orientation, and 
alignment, respectively |2l| . are in Cartesian elements 


1 'il-jk ~ q fiilfijkt 


1'il-jk 9 (fiijfikl SikSji), 


1'il-jk ~ 2 + bikbjl) (39) 

and the corresponding eigenvalues can be expressed by 
6 _j-symbols of angular momentum theory: 


(37) 

(38) 


1 


— 3(2 J e + 1) 


11 k 

Jr Jr J? 


(40) 


Using the decomposition (36), the entire multiple scat¬ 


tering series for the average intensity has been summed 
exactly in Ref. m, including the constraint of transver- 
sality for the far-field photons exchanged between distant 
atoms. In the diffusion approximation, the bulk intensity 
propagator (known as the “diffuson”) for the polarization 
channel e a H > e b reads in Fourier representation 


r ab (q) = £ 


L ab 


k ;—0 


q 2 + {Dtk )- 1 ‘ 


(41) 


Here, the momentum q is the Fourier variable conjugate 
to the position difference, and the mode populations 


f( K ) \ c c -* rp(K-) 

L ab / v tgit-bj t-bktgi-L ji.jk 

ijkl 


( 42 ) 
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depend only on the initial and final polarization. Each 
mode has a diffusive propagator with a characteristic ex¬ 
tinction range l 2 = Dt k . Here, D = l 2 /3r is the (scalar) 
diffusion constant with l the mean free path and r the 
energy transport time m ■ The depolarization times 


A k t 

3(1 -b K X K ) 


(43) 


describe the damping of mode populations. This depo¬ 
larization is due to two effects: First, the transversality 
constraint of free-space propagation between atoms, as 
encoded by the photon propagator eigenvalues 

bo = 1, h = *, b 2 = (44) 


Second, a scrambling of polarization by internal atomic 
Raman transitions, as encoded in the vertex eigenvalues 
A k , Eq. (|40l. 


Only the scalar mode k = 0 that counts the total in¬ 
tensity or number of photons has an infinite lifetime, 


Tq =0, since & 0 A 0 = 1 for all J e , J g . The two other 
modes k = 1,2 have b K X K < 1 for all J g , J e and thus de¬ 
scribe depolarization on the scale of the transport time 


C. Multiple scattering: the correlon 


The previous derivation for the diffuson can be applied 
to the correlation propagator of ( 20 ) (dubbed henceforth 


“correlon”), by substituting the correlation vertex (35) 


for the intensity vertex (33). The correlation vertex de¬ 
composes as 


r( c ) _ \ ' \(c )rp(n) 

1 il;jk — / a k 


K — 0 


il]jk 


(45) 


over the same projectors (37) (39). Bot now there is 


only a single, fully degenerate eigenvalue for all modes 

k = 0 , 1 , 2 : 


\( c ) = M = 2 ^ e + 1 
K es 3(2 J g + 1) ‘ 


(46) 


This is nothing but a complicated way of saying that 
the vertex (351 is proportional to the identity, preserving 


the polarization under Rayleigh scattering. The diffusive 
correlation propagator then takes the same form as (41), 


^ab,a'b' (*0 ~~ <^k a -k b ,k„/-k b , 


t 


(k) 


^o ? 2 + 


(47) 


The prefactor expresses the conservation of total momen¬ 
tum, restored by the ensemble average over atomic po¬ 
sitions, for the four independent momenta entering the 
correlator. The mode occupation factors also depend on 
all four polarizations concerned, 

taia'b ( 48 ) 

ijkl 


In the parallel channels e a i = e a and = e&, one finds 


t^o.b = t'ab as gi ven by (|42j) . The decorrelation times 




r (c) = 

' K 


M PI 


3(1 - b K M eg ) 


(49) 


are non-negative under all circumstances because 0 < 
b K < 1 and 0 < M eg < 1 for all J e , J g because of the 
dipole selection rule | J e — J e \ < 1. 

Only atoms with a non-degenerate ground level (J g = 
0, J e = 1) act as scalar point scatterers, and there 
is no difference between the average intensity and the 
average amplitude squared (as discussed at length in 
the context of coherent backscattering 0), such that 


(T^tab) = (T o 6 )*(t a6 ). Equivalently, the correlation 
vertex and the intensity vertex have the same eigenval- 

(c) - 

ues, A K = = 1. Therefore, from Eq. (25) we recover 

the Rayleigh law of unit relative fluctuations, 


O ab,ab — 1 (Jg — 0)- 


(50) 


But in general, contrary to what is claimed in [9j, we 
find no negative decorrelation times that could signal an 
instability and a build-up of correlations over large dis¬ 
tances. For degenerate dipole scatterers with J g > 0, the 
correlation vertex eigenvalues differ from their intensity 
counterparts, but, as shown by (27), the speckle corre¬ 
lations can only be reduced by the internal degeneracy, 
never enhanced. 


D. Speckle correlations across a slab 


As an application of the preceding formalism, we com¬ 
pute the forward correlations a = af, b = b' measured 
in transmission across a slab of width x = L/l in units 
of the mean free path l, following [9]. The solution of 
the diffusion equation in the slab, for each irreducible 
polarization mode n, is 


r.w = 

7« smh7 re a: 


(51) 


with reduced relaxation rates 


7k = VtJt.K = V^A * 1 - b K ) 


(52) 


for both the intensity (whose eigenvalues X K are given by 
(40)), and the correlation (whose eigenvalues are A ^ = 
M eg ). For conserved modes with = 0, (51) reduces to 
r K (x) = 1 /x. 

Concerning the polarizations, we assume unpolarized 
incident light and no polarization selection at detection, 
which implies that the weights t^ ab = given by (42) 
are averaged to 


J°) = 2 

ab 3 ’ 


i W — o 

b ab ^ U ’ 


J 2 ) _ 

L ab ~ 


1 

3' 


( 53 ) 
























ground-state angular momentum J g 


FIG. 2. (Color online) Decorrelation rate, Eq. (521, as func¬ 
tion of ground-state angular momentum J g (in units of K) 
for the 3 possible dipole transition types | J e — J g | < 1 . The 
hatched area indicates the region 7 q < 0 forbidden by the 
Rayleigh bound. Except for the non-degenerate case J g = 0 
of isotropic dipoles, where the rate vanishes and the Rayleigh 
law applies, the decorrelation rate is always larger than unity, 
signalling a rapid loss of correlations. 


The relative fluctuations, Eq. (25), thus read 


Cab,ab{x) 


/ ^(aQ+lf^A 2 3 

^ 2r 0 Or) + r 2 (:r) ) 


(54) 


In the denominator, and for thicknesses x 1, the dif¬ 
fusive Goldstone mode k = 0 with its algebraic decay 
dominates over the exponentially decaying mode k = 2 . 
In the numerator, bo = 1 > fe 2 = 7/10 implies 7 ^ < 7 ^. 
Thus, the mode n = 0 has the slower decrease and will 
dominate for thick samples. For large optical thickness 
x 1 , we therefore find that the correlations 

C a b,ab(x) ~ x 2 exp{- 27 ^ c) a;}, (J g > 0), (55) 

are exponentially small, with a rate 

7q C) = V^eg 1 - 1 ) (56) 


per mean free path. This decorrelation rate is plotted in 
Fig .[2]as function of J g for the 3 different dipole transition 


types J e ~ Jg = 0, ±1. For J g = J e , one has M eg = 1/3, 
and the prediction 7 ^ = 3y/2 is independent of J g . For 

all J g > 0, the rate 7 ,) per mean free path is larger 
than unity, signalling a rapid loss of correlations. The 
only case where correlations should survive is for isotropic 

(c) 

dipoles, J g = 0, J e = 1, where 7 / = 0, and the Rayleigh 
law (50) applies. 


IV. SUMMARY AND OUTLOOK 


We have calculated analytically the short-range speckle 
correlations measured by light scattering across a dilute 
cloud of atoms with internal degeneracies, as proposed 
by 0. These correlations are due to the interference 
between certain pairs of disjoint path amplitudes. This 
interference can only occur in the absence of Raman scat¬ 
tering events, which act as path markers. In agreement 
with our initial claim tun. we find that these correlations 
relative to the average intensity are always smaller than 
unity and cannot exceed the Rayleigh limit, which is only 
recovered for atoms with a non-degenerate ground state 
(J g = 0). In transmission across an optically thick slab, 
the correlations are even predicted to be exponentially 
small. Of course, the exponential reduction (55), within 
the strict protocol proposed in Ref. 0 , does not preclude 
the enhancement of spectral fluctuations by internal de¬ 
generacies, as predicted for an entirely different protocol 
in Ref. [22]. 


The internal degeneracy of atomic point scatter¬ 
ed is thus found to reduce generic mesoscopic effects 
such as diffusive intensity correlations and coherent 
backscattering. As a possible extension of this work, it 
would be interesting to analyze whether more long-range 
correlations—due to paired paths with scattered in com¬ 
mon [ 2 ,— are affected in a similar way. 
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